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ABSTRACT 

The size distribution and total mass of objects in the Oort Cloud have important implications to the 
theory of planets formation, including the properties of, and the processes taking place in the early 
solar system. We discuss the potential of space missions like Kepler and CoRoT, designed to dis- 
cover transiting exo-planets, to detect Oort Cloud, Kuiper Belt and main belt objects by occultations 
of background stars. Relying on published dynamical estimates of the content of the Oort Cloud, 
we find that Kepler s main program is expected to detect between and ~ 100 occultation events 
by deca-kilometer-sized Oort Cloud objects. The occultations rate depends on the mass of the Oort 
cloud, the distance to its "inner edge", and the size distribution of its objects. In contrast, Kepler is 
unlikely to find occultations by Kuiper Belt or main belt asteroids, mainly due to the fact that it is 
observing a high ecliptic latitude field. Occultations by Solar System objects will appear as a pho- 
tometric deviation in a single measurement, implying that the information regarding the time scale 
and light-curve shape of each event is lost. We present statistical methods that have the potential 
to verify the authenticity of occultation events by Solar System objects, to estimate the distance to 
the occulting population, and to constrain their size distribution. Our results are useful for planning 
of future space-based exo-planet searches in a way that will maximize the probability of detecting 
solar system objects, without hampering the main science goals. 

Subject headings: Oort Cloud — Kuiper Belt — comets: general — techniques: photometric — 
stars: variables: other 



1. INTRODUCTION 

Oort (1950) postulated the existence of a cloud of 
comets orbiting the Sun with typical semi-major axes, 
fl, of ~ lO'* AU. This cloud is required to explain the ex- 
istence of long period comets and the fact that a con- 
siderable fraction of these comets have orbital ener- 
gies concentrated in a narrow range, corresponding to 
1/a « 10^^ AU^^. However, to date no direct observa- 
tion of objects in the Oort Cloud exist"*. 

Dynamical simulations suggest that the Oort Cloud 
was formed by the ejection of icy planetesimals from the 
Jupiter-Neptune region by planetary perturbations (e.g., 
Duncan et al 1987; Dones et al. 2004). On orbital time 
scales, galactic tides and passage of nearby stars raised 
the perihelia of these comets above the region of the 
giant planets influence. Comets with semi-major axes 
above - 2 X lO'^AU are loosely bound to the Sun, and 
tides by the galactic potential (e.g., Heisler & Tremaine 
1986) as well as impulse by passing nearby stars, may 
send these comets to the inner solar system where they 
reveal themselves as long period comets. In contrast, 
Oort Cloud objects with a < 2 x 10"^ AU (' 'inner Oort 
Cloud"; e.g.. Hills 1981, Bailey 1983) have more sta- 
ble orbits. We note that the inner "boundary" of the 
Oort Cloud, rmin, is set by perturbations due to the gi- 
ant planets and the primordial stellar neighborhood of 
the Sun. This rmin is estimated to be in the range of 
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1,000AU to 3,000 AU (Duncan et al. 1987; Dones et 
al. 2004; Fernandez 1997). Below 

''min the number of 
comets steeply falls due to planetary influence. Dynam- 
ical simulations suggest that the spatial density of Oort 
Cloud objects in the 3,000-50,000 AU range falls as r", 
with a = -3.5 ± 0.5 (e.g., Duncan et al. 1987). 

The current flux of long period comets calibrated by 
dynamical simulations suggest that the outer Oort Cloud 
(i.e., fl > 2 X 10^ AU) contains ^ lO'^ comets with nu- 
clei absolute planetary magnitude^ brighter than 16 mag 
(e.g., Heisler 1990; Weissman 1996). This magnitude 
roughly corresponds to objects with radius of 2 km (as- 
suming 4% geometric albedo). The dynamically more 
stable inner Oort Cloud contains about 2-100 more 
comets than its outer counterpart (e.g.. Hills 1981; Dun- 
can et al. 1987; Fernandez & Brunini 2000; Dones et al. 
2004; Brasser et al. 2008). 

Detecting Oort Cloud objects, and estimating their 
size distribution, will improve our understanding of the 
mass of the solar-system planetary accretion disk, and 
the dynamical processes in the young solar system. 
Moreover, it may bear some clues to the stellar density 
at the region in which the Sun was born (e.g., Fernandez 
1997; Fernandez & Brunini 2000; Brasser et al. 2006). 
However, these objects are beyond the direct reach of 
even the currently largest planned telescopes, and some 
indirect methods (e.g., astrometric microlensing; Gaudi 
& Bloom 2005). 

Bailey (1976) suggested that small Kuiper Belt ob- 
jects (KBOs) as well as Oort Cloud objects could be 
detected by occultations of background stars. Indeed 
several surveys are looking for KBO occultations (e.g.. 

Defined as the magnitude of an object observed at opposition and 
at 1 AU from the Sun and Earth. 
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Roques et al. 2006; Lehner et al. 2009; Chang et al. 
2006, 2007; Bianco et al. 2009; SchUchting et al. 2009). 
However, to date there is only one reported occupation 
of a Kuiper Belt object (Schlichting et al. 2009). 

In this letter we discuss the potential of space missions 
designed to search for transiting exo-planets (e.g., Ke- 
pler, CoRoT) to detect occupations of background stars 
by small Oort Cloud and Kuiper Belt objects. 



2. STELLAR OCCULTATIONS BY SMALL SOLAR SYSTEM 
OBJECTS 

Space telescope designed to look for exo-planet tran- 
sits usually have exposure times, 4xp. longer by orders of 
magnitude than the typical duration of an occupation of 
a background star by a small solar system object, fjur- 
The relative flux decrement of such an occupation A, 
integrated over the exposure time, is diluted by a fac- 
tor oc fdur/fexp- However, these satellites have superb 
photometric accuracy. Their noise level per exposure is 

(7 = Sts:^J'^ ~ lO""* — 10~^, where 5 is a normalization 
parameter^. Therefore, we may be able to detect these 
minute occupations. Here, we adopt a detection crite- 
rion for an occupation of A > na, where n is the signal- 
to-noise ratio. Throughout the paper we adopt n = 8. 

For an occultation by objects at a given distance from 
Earth, r, the maximal value of A (obtained for an optimal 
impact parameter^ b) increases monotonically with the 
object size. We define 7?min as the radius of the small- 
est object that can be detected, and correspondingly 
6min = Rmin/ f. Furthermore, we can define the sum of 
impact parameters for which objects with a radius R are 
detectable, bz{R,na) =2 J^&[A{R,b) - na]db, where 
© is Heaviside step function. According to these defini- 
tions bz{R < /?min,«cr) = 0. 

The size distribution of KBOs and Oort Cloud ob- 
jects per unit solid angle is parametrized by a ^o^er- 
law, i/A'(/?)/i//? = A'> 1 km ( 1 km) - 1 - 1 ) 1 km) or 
by a broken power-law. Here, A^>ikm is the total num- 
ber of objects with R> \ km per unit area, and q is the 
size distribution power-law index. For Oort Cloud ob- 
jects AA^>ikm ~ lO'^, where A is the surface area of the 
celestial sphere. For KBOs larger than /?brcak ~ 45 km 
q = 4.5 (Bernstein et al. 2004; Fuentes & Holman 2008; 
Fraser et al. 2008), and (7 w 3 — 4 for smaller KBOs (e.g, 
FarineUa & Davis 1996; Pan & Sari 2005; Schlichting et 
al. 2009). For Oort Cloud objects q is unconstrained by 
observations (see however Goldreich et al. 2004). 

The number of events in a survey that ob- 
serves A'* stars for a duration T is then A'ev — 
A^*T Jg^^iJ.bz{R,na)j^dR, where n is the angular ve- 
locity of the occulting objects on the sky as seen from 
Earth. Here all objects are at the same distance r, all 
stars have the same magnitude and angular radii and n is 
constant. When this is not the case, it is straightforward 
to integrate over the distribution of these properties. 

The functional forms of tdm, A, /?min and thus bz{R) 
depend on the ratios between the three angular scales 
in the problem. These are the object angular radius. 

For Kepler (CoRoT) S = 0.0022 (0.045) for 1 s integration of a 
12th magnitude star. Tliis sensitivity parameter depends on the star 
magnitude m as iQ+O-^im-n) ^ 

^ In the case of diffraction occultation it is not at = 0. 



©obj =R/r;tiie stellar angular radius^ 0*; and the angu- 
lar Fresnel scale^ 6p=F/r, with F = y^Xr/2, and A the 
wavelength of observation. The duration of the eclipse 
is tdui ~ 20max/Al, where 0max = max(0obj, 0F, 6*). Ta- 
ble 1 gives the functional form of T^min in the different 
"asymptotic" cases (denoted by A-E), where the param- 
eters are normalized to those of a typical Oort Cloud 
object and the Kepler sensitivity. 

The cases where 0obj ^ 0*,dp (case A) and 
0* dobj,dp (cases B-C) represent geometric oc- 
cultations. In these cases we approximate bz{R) = 
max{0obj,0*}(ll-7 -q)/13 for any R > R^in- Here, 
the (^-dependent factor (i.e., [11.7 — (j'J/lS) is aresult of 
the assumed circular shape of the object. We obtained 
this approximate correction factor by numerically inte- 
grating geometrical occupations by circular shaped ob- 
jects^". In the diffractive regime (Op ^ 0*,0obj; Cases 
D-E) we numerically calculated diffractive light curves 
using Eq. Bl in Roques et al. (1987). We find that bz 
can be well approximated by 



bz{R,no) S 



max{0obj,0*} 



exp 



^ no t, 
4 ?du 



9p_ 
Oobj 



(1) 

The scaling of Eq. 1 can be understood as follows. The 

occultation "shadow" pattern of an object can be de- 
scribed as a set of bright and "dark" concentric rings 
with the central ring (around the "Poisson peak") be- 
ing dark with an angular radius of dp. The area in aU 
the rings is constant (~ 0p) so the angular width of the 
rings is inversely proportional to its angular distance 
from the center (e.g.. Fig. 1 in Nihei et al 2007). The 
'depth' of the decrement in the flux in the dark rings 
is PS (^obj/^j)^ so the value of A for the central ring 
occultation is Ac ^ (tdm /tcxp) (Oobj /dp)^. This decre- 
ment remains similar in all the 'dark' rings as long as 
the ring width is smaller than max{0*, 0obj}- Therefore 
bz{R,nG<Ac) - 9^/ max { , } . 

For each of the five cases we use the approximations 
described above and give formula (Table 2) for the ex- 
pected number of events as a function of the parameters 
n, S, texp, d*, r, A and Vrei (= ^lr). Here Vrei is the rel- 
ative velocity between the occulting object and Earth, 
projected on the plane perpendicular to the Earth-object 
direction. Since the object size distribution is steep, the 
integration over R is dominated by the smallest objects, 
and is carried to infinity (rather than the radius at which 
the condition defining the "case" is violated). 

3. OCCULTATION RATES FOR kepler AND COrot 

To estimate the observed rate of events we need to take 
into account the radial distance distribution (discussed in 
§1), the relative velocity (vrei) probabiUty distribution, 
the ecliptic latitude distribution, and the stellar angular 
size and magnitude co-distributions. 

3.1. Stellar angular size and magnitude distribution 

9t = Te^^lO""-^*'''"', where Tc and are the effective tempera- 
ture and bolometric magnitude of the star, respectively. 

9 For r = 3000 AU and A, = 5000 A, F=10.6 km and = 4.9 x lO'^ 
arcsec. 

This integral do not have an analytical solution. 



Table 1. 

JECTS 



Kepler and the Oort cloud 

OCCULTATION REGIMES AND THE MINIMUM RADIUS OF DETECTABLE OB- 
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Case Condition 



[km] 



A 


6obj S> 9» , Op 


7.9 


B 


e. > Qobj > % 


16.6 


C 


e. > Qf > Sobj 


14.0 


D 


Of > Sobj > e. 


4.6 


E 


Of > e. > Sobj 


4.6 



K 5 

8 0.0022 



'cxp 



,1/2 



1/2, 



1/2 



900s; aokms-' 

8 0.0022 j \<Mh) l,30km,s-l ) l,3000ALlj 1,1.6x10-5" ) 

8O.OO22J l^SOOs^ l,30kiiis-l / l^SOOOAU^ I, 1.6x10-"' ) 

- \l/2/, n1/4/ „ \l/2/ \I/4 , 

8 0.0022 , 



I, 900s } 



aOkms 



, 3000 AU , 



' n S 

, S 0.0022 ; I 900 



1/2/ , y/4 



5000 A 



1/4 



V 5000 a/ 



Note. — The normalization of the star angular radius, 6,, corresponds to S = 35 km at r = 3000 AU. The difference between case B and case C 
is that in case C we replace 9„\,j by 2'/^9oi,j. This small coiTection is required in order to take into account the increased cross section of diffractive 
occultations (estimated numerically). All the formulas assume that fexp > !dur- 



TABLE 2. Formulas for the expected number of occultation events in the 

DIFFERENT ASYMPTOTIC REGIMES 



Case Wev 



A 
B 
D 
E 



en inO''0(3.5-g) 11.7-g 9-1 AjV^ltni iV. 1: { n S \^ ( 'exp ^ 
13 q-2 10" ToS" 3yr 1, 8 0.0022 j \9ms j 



„ inl.22(3.5-al 11.7-,; AiV>itm N, x f n S ^ 
^■^ 13 10" 105 3yr I 8 0.0022 . 



^sn V 1 n0-66(3.5-(,) 3.5(,7-l) AjV^i^n, N, 



i n S ] 
3yr I 8 0.0022 J 



V 30kms- 

' ia_ \ ( v,d 

,900sj l,30kms- 

, 900s ) 



i-1 



. 1.6x10-5" , 



, 30kms 



V 5000 A/ 



36x100.66(3.5-,)!^ 



N, r I II S 



10" 3yr I 8 0.0022 



(solMl) (safe) (30km's-') (3OOOAU) (l.6x'l'o-5") (^^) 



Note. — Case C is like case B but multiplied by (1/2)'' A is the area of the celestial sphere. Note that these formulas are not continuous 
and they provide a better approximation in the "asymptotic" cases. All the formulas assume that t^^p > 'dur- 



In Figure 1 we present the estimated angular radii dis- 
tribution of target stars observed by Kepler and CoRoT. 
To estimate this we queried the guide star catalog (GSC; 
version 2.2; Lasker et al. 2008) for all stars within 2 deg 
from the center of the Kepler (CoRoT) field, that have 
7?-band magnitudes between 9 and 14 (12 and 16), and 
B R colors between 0.59 and 3.4 mag. This color range 
roughly corresponds to the spectral range of Kepler tar- 
gets. We converted the color index of each star to an 
effective temperature and a bolometric correction by fit- 
ting the color indices with those obtained from synthetic 
photometry of black-body spectra with different tem- 
peratures (e.g., Ofek 2008). The synthetic photometry 
(Poznanski et al. 2002) was preformed using the trans- 
mission curves of the filters used in the GSC (Moro & 
Munari 2000). 

3.2. Relative velocity distribution 

Given the unit vector in the direction of a target f, 
the speed of the observer relative to the occulting ob- 
ject V, and the heliocentric velocity unit vector of the 
observer' ' v, the approximate relative velocity, of an ob- 
server in an Earth trailing orbit, projected on the plane 
perpendicular to f is 



Vrei = vy 1 - (y • r)2 ^ v© y 1 - cos2 (A - Ae) cos2 j3 , 

(2) 

" Object velocity is neglected. 




Fig. 1. — Approximate stellar angular radius distributions for Ke- 
pler (black line) and CoRoT (gray line). The vertical red dashed 
lines represent the Fresnel radii at selected distances (see labels) where 
A = 5000 A. 



where A and )3 are the ecliptic longitude and latitude 
of the observed stai". A© is the ecliptic longitude (as 
measured from the Sun) of the observer, and v© « 
29.8 km s^ ' is the heliocentric speed of the observer Fi- 
nally, assuming circular orbit (i.e., the probability distri- 
bution of A© is uniform) the probability distribution of 



3.3. Ecliptic latitude distribution 

The Oort Cloud is predicted to have a roughly spheri- 
cal shape. However, the Kuiper Belt and the asteroids 
main belt are concentrated toward the ecliptic plane. 
Using the inclination distribution of KBOs estimated 
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Table 3. Approximate occultation rates 



Parameters Rates (T = 3 yr) 



Kepler CoRoT 

A'>lkm '■min q /j, = 6 X l(r* /g = 1 

[deg-2] [AU] = 10^ J3 « 66° A?* = 10"*, jS =i 0° 

2.4x 10** 3000 3.0 30 0.1 

2.4 X 10** 3000 3.5 8 4 x 10"^ 

2.4 X 10* 3000 4.5 0.7 1 x lO"' 

2.4 X 10* 1000 3.5 60 0.03 

2.4 X 10* 5000 3.5 3 2 x lO'^ 

1.1 X lOVe 40 3lF 2 X 10"* 3 X 10"^ 

7.3xl0Ve 40 3.5 5x 10"'' 4x10"^ 

4.9 X 10^ /e 40 4.0 5 x IQ-^ 0.02 

47/e 4 2.3" 1 X 10"'* 0.2 



mates in Table 3 are an approximation to the actual rate. 
This is mainly because our formulas represent asymp- 
totic approximations, while in some cases two out of the 
three angular scales, Oq^, 9p, 0*, maybe be similar. 

4. VALIDATION AND DEGENERACY REMOVAL 

Contrary to transits and grazing transits by extra-solar 
planets, solar system occultations will appear in a single 
photometric measurement and they happen only once. 
Therefore, verification of each event by itself will re- 
quire a thorough understanding of the noise properties 
of the photometer. Kepler will have about 1 .5 x lO'*^ 15- 
min independent photometric measurements. Assuming 
purely Gaussian noise, the 8- (7 detection threshold we 
adopted corresponds to a probability of ^ 10"'^. How- 
ever, even small deviations from a Gaussian noise may 
impair our confidence in the astrophysical nature of the 
events. 

Here we present two tests that can be used to check if 
the occultations signature are due to real events or due 
to some sort of uncharacterized noise. In addition, these 
tests can be used to statistically measure q and to distin- 
guish between Oort Cloud objects and KBOs. 

The occultation rate depends on the relative veloc- 
ity v'rci, with a somewhat different dependency in each 
regime. Therefore, if a large number of occultations are 
detected, a simple sanity test is to look for the depen- 
dency of the events rate on Vrei. We note that since Ke- 
pler is observing at high ecliptic latitude the Vrei modula- 
tion will be small (Eq. 2) compared to the maximal mod- 
ulation which will be seen by an instrument observing 
the ecliptic. Moreover, in all the cases (A-E), the func- 
tion A'ev depends on q. Therefore, the amplitude of the 
modulation is related to q and the occultation type (i.e., 
A-E). Knowing the occultation regime (see below), we 
can measure q (i.e., removing the degeneracy between q 
andA/ov). 

Furthermore, as shown in Table 2 the occultations rate 
depends on S, which depends on the stellar magnitude 
and on 0*, which in turn is a function of the star mag- 
nitude and effective temperature. We expect bright and 
angularly small stars to have higher chance of produc- 
ing detectable occultations. Moreover, this dependency 
is different for the various occultation cases. This can 
be used to further test the nature of the occultations, 
by comparing the distribution of magnitudes and effec- 
tive temperatures of the stars showing occultations to 
those of the entire observed stellar population. Most 
importantly, for a Kepler-like observatory, KBO occul- 
tations will be mostly in the diffractive regime, while 
Oort Cloud object occultations will be in the geomet- 
ric regime (Fig. 1). Therefore, the magnitude and 6^, of 
occulted stars by KBOs and Oort Cloud objects will fol- 
low different co-distributions. Therefore, we can statis- 
tically distinguish between Oort Cloud and KBO occul- 
tation events. However, quantifying these effects and es- 
timating the number of occultations needed to carry out 
the suggested analysis requires a more realistic estimate 
of the events rate (i.e., using full numerical integration 
of light curves) and a full exploration of the parameters 
space. 

To summarize, we show that Kepler and similar space 
missions can detect deca-kilometer objects in the Oort 
Cloud. We present statistical methods to verify that the 



NoiK. — The first block refers to Oort Cloud objects assuming 
(•max = 50000 AU, a = -3.5. The second block refers to KBOs as- 
suming ("max = 40 AU. The third block refers to main belt asteroids 
assuming rn,ax = 4AU. The cumulative surface density of KBOs near 
the break radius, ^break ~ 45 km, is 5.4 deg"^ at the ecliptic (Fuentes 
et al. 2009). Therefore, at the ecliptic we adopt for KBOs, A'>ikm of 
1.1 X 10", 7.3 X 10*, and 4.9 x 10^ deg^^, for q =3, 3.5, and 4, re- 
spectively. We note that these values are consistent with the findings 
of Schlichting et al. (2009). A'>ikm and q for the main belt are taken 
from Ivezic et al. (2001). We set the ecliptic latitude surface density 
correction /j to 6 x 10"' and 1 for ji f» 66° and ji ~0, respectively, 
for both KBOs and main belt asteroids. In all the cases we assume 
eight sigma detection threshold (i.e., n = 8). 

^For KBOs we assume a broken power-law with q = 4.5 above 
R=45km. The actual rate calculation is preformed by normalizing 
the Equations in Table 2 to 45 km, and using different q for above 
and below 45 km. 

''For main belt asteroids we assume a broken power-law with q = 4 
above R=2.5 km, and q = 2.3 below this radius. 



by Elliot et al. (2005) the expectation probability to 
find a KBO at an ecliptic latitude of ^ 66° is about 
fe^6x 10"^ times smaller than that at the ecliptic. 
This relative probability was estimated by convolving 
the inclination probability distribution of KBOs with the 
echptic latitude distribution for a given inclination (i.e., 
ElUot et al. 2005 eqs. 40 and 34, respectively). We note, 
however, that this number is highly uncertain. 

3.4. Rate estimates 

To estimate the rate of occultation events we integrate 
the appropriate formulas in Table 2 over Vy^u r, 0« and 5 
co-distributions. For each combination of parameters we 
integrate over the most appropriate occultation channel 
(i.e., cases A-E). In Table 3 we present the estimated 
rates predicted for Kepler and CoRoT for the Oort Cloud 
(upper block), Kuiper Belt (middle block), and main belt 
(lower block) object occultations. The predicted rate for 
Oort Cloud objects in Table 3 in case of q = 3.5 and 
rmin = 3000 AU is lower by about an order of magnitude 
compared to the pre-factors in Table 2. This is mostly 
because Table 2 is normaUzed for all objects being at r = 
3000 AU while in Table 3 the same number of objects 
is spread between rmin and r^ax, with a spatial density 
distribution proportional to r~^-^. 

These rate estimates suggest that if our current 'best 
guess' regarding the properties of the Oort Cloud are 
correct, then it may be detectable by Kepler. Regardless, 
Kepler will be able to put the best constraint (so far) on 
the content of the Oort Cloud. We note that the rate esti- 
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occultations are real rather than due to uncharacterized 
noise. In addition we suggest that it will be possible to 
statistically measure their size distribution {q), and the 
dominant occultation regime, and therefore differentiate 
between Oort cloud objects and KBOs. 

Considering Kepler capabilities and reasonable Oort 
Cloud parameters, we find that Kepler may detect to 
~ 10^ stellar occultations by Oort Cloud objects. We 
find that Kepler is unlikely to detect Kuiper Belt objects 
(mainly because of its high ecliptic latitude pointing). 
Moreover, CoRoT is unlikely to detect Kuiper Belt or 
Oort Cloud objects. However, the exact properties of 
the Oort Cloud and the Kuiper Belt are not well known. 
Therefore, such searches are warranted. 



The analytical treatment we present in this paper is 
useful in order to maximize the efficiency of existing 
(and future, e.g., PLATO) experiments to detect trans 
Neptunian objects. However, this treatment is accurate 
only in the asymptotic cases, and numerical calculations 
are needed in order to give more precise predictions. 
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